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Abstract 

We present a theoretical study on the band structures of the electron constrained to move along 
triply-periodic minimal surfaces. Three well known surfaces connected via Bonnet transformations, 
namely P-, D-, and G-surfaces, are considered. The six-dimensional algebra of the Bonnet trans- 
formations [C. Oguey and J.-F. Sadoc, J. Phys. I France 3, 839 (1993)] is used to prove that the 
eigenstates for these surfaces are interrelated at a set of special points in the Brillouin zones. The 
global connectivity of the band structures is, however, different due to the topological differences of 
the surfaces. A numerical investigation of the band structures as well as a detailed analysis on their 
symmetry properties is presented. It is shown that the presence of nodal lines are closely related 
to the symmetry properties. The present study will provide a basis for understanding further the 
connection between the topology and the band structures. 

PACS numbers: 73.22.-f, 61.46. +w, 68.35.Bs 
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I. INTRODUCTION 



For nearly two decades intense efforts have been devoted to the studies of periodic curved 
surfaces as topological objects to describe certain structured materials. Several conceptu- 
ally different classes, such as zero-potential surfaces, nodal surfaces, and periodic minimal 
surfaces, often share similar properties concerning topology and symmetry.— ^ Such a surface 
can be constructed as a space partitioner which is invariant under space group operations, 
hence the space group symmetries are the key ingredients of them. In reality, on the other 
hand, such surfaces are realized as nanometer to micron scale interface structures in ly- 
otropic liquid crystals or cell membranes. Furthermore, recent progress in the synthesis of 
nano-structured materials prompts us to consider the possibility of using these materials as 
the templates for electronic devices embracing interesting geometrical properties. 

The present paper deals with a fundamental problem on the electronic structures of 
materials characterized by periodic curved surfaces. The theoretical basis of the quan- 
tum mechanics of electrons confined within thin curved films has been developed by several 
authors . 3 i 4 i 5 i 6 i 7 i 8 They showed that when the electron is strongly constrained onto a smoothly 
curved surface, the quantization of the motion perpendicular to the surface results in an ef- 
fective potential energy depending on the local curvatures along the surface. Since this 
potential energy is attractive, the electron can be weakly bounded around a curved region 
of the surface (geometrically induced bound states). Theoretical interests also arise in re- 
gard to quantum chaos, since the constrained motion on a surface with non-zero intrinsic 
curvatures entails classical chaos. 9 These studies, however, deal primarily with the local ef- 
fects of curvatures. The roles played by the symmetry and topology of the surface on the 
determination of electronic structures have largely remained intact, apart from that Aoki et 
al. numerically studied the band structures of triply periodic minimal surfaces recentl y 10 ! 11 
The aim of the present paper is to perform a more thorough investigation into the electronic 
structures of similar periodic curved surfaces. 

In SecJnJ we provide a mathematical background for periodic minimal surfaces and in- 
troduce the three basic surfaces in the same Bonnet family, i.e., P-, D-, and G-surfaces, 
which will be studied further. In Se dllll the translational subgroups of these three surfaces 
are described as orthogonal projections of an identical six-dimensional (6D) translational 
subgroup onto 3D subspaces. A numerical investigation of the band structures of the three 
surfaces is presented in Sec llVI following a brief summary of the basic formulation of the 
one-electron problem constrained onto a surface embedded in 3D space. In SecJVjthe alge- 
braic information given in Sec JIIII is used to prove certain interrelations between the band 
structures of the three surfaces. More precisely, we propose the idea of "band intersections" 
to show that the energy eigenvalues at a certain set of wave vectors coincide among the 
three surfaces. In Sec JVIl the symmetry properties of energy bands, i.e., the corresponding 
irreducible representations (IR's) of the space groups, are analyzed. Based on the IR's, 
several features of the band structures and eigenstates are discussed, including the band 
sticking phenomena and the occurrence of nodal lines. In particular, a general connection 
is established between the symmetry properties (or IR's) of the eigenstates and the types 
of nodal lines entailed by symmetry. The exhaustive lists of IR's with such nodal lines are 
shown. Further discussions are given in the last section. 
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II. PERIODIC MINIMAL SURFACES 



Curved surfaces are naturally realized as interfaces, and their formation is dominated 
by surface energies (e.g. surface tension). Like soap films on closed wire frames, they 
often resemble minimal surfaces, that is, surfaces with vanishing mean-curvature. Minimal 
surfaces are privileged mathematically since they have an analytic representation called the 
Weierstrass-Enneper (WE) representation : 12 ! 13 a minimal surface is represented by the real 
part of a complex vector z in a 3D complex space C 3 (~ E e ) multiplied by a phase factor: 

r(w) = Re (e ia z(w)) , (1) 



z(w 



J F(C)(1 - C 2 )dC, J «F(C)(1 + C 2 )dC, J 2F(C)CdCj , (2) 

where I is an overall scale factor, a the Bonnet angle, w = q 1 + iq 2 , and F(() a complex 
analytic (or holomorphic) function associated with the Bonnet family of minimal surfaces.— 
Hence the surface is represented as a multiple- valued function of the complex variable w, 
where different integration paths give distinct points on the surface. The two real variables, 
q 1 and q 2 , serve as the isothermal parameterization of the surface. 

Minimal surfaces with periodicity in three independent directions are called triply periodic 
minimal surfaces (TPMS's). The importance of TPMS's as topological objects describing 
interface materials like self-assembled lyotropic liquid crystals has been recognized for nearly 
two decades.— i 15 i 16 Compared to other classes of periodic surfaces like zero-potential surfaces, 1 
TPMS's are especially advantageous for physical applications because Eq.Q provides a 
parameterization of the surface for incorporating it into written equations. In the following, 
we focus on the classical examples of Schwarz's P- and D-surfaces^ and Schon's G-surface 
(gyroid)^ 8 - which belong to the cubic system. The corresponding WE function, F((), is given 
by 

^(0 = 1 (3) 

y/1 + 14C 4 + C 8 

with eight branch points of order 2 at £ = (\/3± 1) / y^e 2 ™/ 4 (n = 0, 1, 2, 3) on its Riemann 
surface.— For each of the three surfaces, the Bonnet angle a, has to take a proper value as 
will be given below. Multiple surfaces in the same Bonnet family are intrinsically identical, 
or isometric, apart from the global topology. The surfaces are obtained from one another 
via Bonnet transformation, namely the change of a. 

We may write z(w) = u(w) + iv(w) with u(u>) and v(w) being two real vectors, hence 
r(w) = u(w)cosa — v(w)sina. Alternatively, one may introduce a 6D vector, H(iy) = 
(u(w), v(w)), and the minimal surface r(w) is given as the orthogonal projection of S(u>) 
onto a 3D subspace specified by the hyper- rotation angle a. P- and D-surfaces are given 
by the Bonnet angles, ap = and old = 90°, hence they are given by u(w) and — v(w), 
respectively. Accordingly the 6D space E 6 embedding E(w) is divided into two 3D subspaces 
as Eq = Ep © Ed, where Ep and Ed are the embedding subspaces for P- and D-surfaces, 
respectively. 

It is convenient to define an orthogonal transformation of E e by the 6D orthogonal matrix 

„, \ / I cos a — I sin a \ ,. s 
0(a) '.— I t • -p J (4) 
v ' \ Ism a I cos a / v ' 
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with I being the 3D identity matrix; the two 3D subspaces, Ep and Ed, are mixed by such 
a transformation. Then Eq.Q can be written concisely as r(u>) = ITpO(a)S(u>), where lip 
is the orthogonal projector onto Ep defined by (I O) with O being the 3D zero matrix. 
G-surface is obtained by choosing the Bonnet angle as «g = arctan(r/s) ~ 51.985°. Here 
the constants r and s are defined by r — is = 2e~ ni/6 K(k 2 ), in which K(k 2 ) stands for 
the complete elliptic integral of the first kind and k = e~ m ^ 3 . Their numerical values are 
r m 2.156 and s ~ 1.686. Hence, G-surface is just the orthogonal projection of E(u>) onto 
a 3D subspace Eq which is inclined in the 6D space by the hyper-rotation angle olq. The 
relevant orthogonal projector is given by = J\pO{olg\ One may also find that the 
orthogonal projector for D-surface is given by = HpO{oid) = — (O I). 

III. CRYSTALLOGRAPHY OF BONNET TRANSFORMATIONS 

With each of the three TPMS's a pair of space groups, G/H, is associated; G includes 
the symmetry elements to reverse the sides of the surface, while H does not. The group H 
is a subgroup of G. The pairs for P-, D-, and G-surfaces are Im3m/Pm3m, Pn3m/ Fd3m, 
and Ia3d/ I4i32, respectively. It is important to note that the Bravais lattices for P- and 
D-surfaces depend on whether or not to identify the sides of the surfaces. The electronic 
Hamiltonians for the three TPMS's to be studied have the symmetry of the group G. 

The WE formulae for the three surfaces reveal that the fundamental patches of the 
surfaces, corresponding to the asymmetric unit of the relevant space groups G, consist of 
(2,4,6) triangles . 20 ' 21 More precisely, the fundamental patch of P-surface as well as that of 
D-surface coincides with a (2,4,6) triangle. That of G-surface, on the other hand, consists of 
two triangles and is twice as large as those of the other two (c.f. FigJTJ). The two triangles 
for the case of G-surface are not congruent with each other in the embedding space, so that 
one of them cannot be reproduced via any space group operation from the other. However, 
they are congruent intrinsically, which implies that G-surface has an additional "intrinsic 
symmetry" which cannot be reproduced by the space group. The surfaces are thus tessellated 
into (2,4,6) triangles. 

A close relationship between the three TPMS's and the (2,4,6)-tessellation of the hy- 
perbolic plane, i.e., the non- Euclid plane with a constant Gaussian curvature, was pointed 
out by Sadoc and Charvolin.— They proved that the translational subgroup of the (2,4,6)- 
tessellation of the hyperbolic plane, where the unit cell is identified as a dodecagonal region 
consisting of 96 triangles as shown in Fig|21 derives the translational subgroups of the three 
TPMS's. The corresponding Bravais lattices are simple cubic (s.c), face centered cubic 
(f.c.c), and body centered cubic (b.c.c.) for P-, D-, and G-surfaces, respectively. The unit 
cells of the TPMS's can be taken as dodecagonal patches embedded in three dimensions, as 
shown in Figd corresponding to the dodecagonal region of the hyperbolic plane. Note that 
these unit cells are only the primitive unit cells for the space groups, H . If the two sides of a 
surface are indistinguishable, the space group will be G whose primitive unit cell is reduced 
by half for P- and D-surfaces, and the Bravais lattices become b.c.c. and s.c, respectively. 

The 6D crystallographic description of S(u>) was thoroughly investigated by Oguey and 
Sadoc,— who established that it is a periodic structure. The unit cell of the relevant 6D 
lattice A is associated with the aforementioned dodecagonal patch lifted into the 6D space. 
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FIG. 1: The unit cells of the three TPMS's obtained by projecting the 6D unit cell (corresponding 
to the dodecagonal region of the hyperbolic plane) onto the relevant 3D subspaces Ep, Ep>, and 
Eg- Bounding boxes are shown to indicate the relative orientation to the cubic unit cell. The 
fundamental patch consists of either a single triangle denoted as A (for P- and D-surfaces) or two 
triangles denoted as A and A' (for G-surface). Representative high-symmetry geodesies of type a, 
f3, and 7 are also indicated on the surfaces. For P-surface these correspond to the two-fold axis U yz 
(a), the mirror planes a xy (/3) and a z (7). For D-surface, these are the mirror plane a xy (a), the 
two-fold axes U zx (/?) and U y (7), respectively. For G-surface, no symmetry element corresponds 
to these geodesies. 

The primitive lattice vectors of A are given by the column vectors o& 

/ r r — r \ 

A = (fx, t 2 , t 3 , f x , f 2 , T 3 ) = i 

\ s -s U s s U / 

The translational symmetry of each of the three surfaces is given as the 3D lattice, Ax, with 
X = P, D, or G. Then it is given as Ax = IlxA. Specifically, Ap (or A^) is generated by 
the six column vectors of the upper (or lower) half of A, so that it is a s.c. (or f.c.c.) lattice 
with lattice constant ap := £r (or an '■= 2£s). Similarly, Aq is generated by the six column 
vectors of the matrix 

/ 1 -1 1 -2\ 
n G A = £6 1 1 -1 -2 (6) 
\-l 1 1 0-20/ 

with b = rsj \Jr 2 + s 2 , so that Aq is a b.c.c. lattice with lattice constant ac ■= 2£b. 
In each case, the unit cell is given by the projection of the dodecagonal patch onto the 
relevant subspace (see FigHJ). The relationship between A and Ax (X = P, D, and G) is 
schematically represented in FigEl 

Let us consider the reciprocal lattice A* of A. To do so, we shall take the transposed 



r — r r 

r r — r 

s — s s s 

— s s s s 



(5) 
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FIG. 2: The dodecagonal region of the hyperbolic plane, shown within the Poincare disc, cor- 
responds to the dodecagonal unit patches of P-, D-, and G-surfaces. The identification of the 
opposite edges yields a genus-3 torus corresponding to the unit cell of E(w) with opposite sides be- 
ing identified. The tessellation into (2, 4, 6) -triangles corresponds to the division of the TPMS's into 
fundamental patches. There are three types of geodesies passing through high-symmetry points: 
a, passing through six- and four- fold symmetry centers alternately, /3, six- and two- fold centers, 
and 7 four- and two-fold centers, respectively. Only a few representatives of each type is indicated. 
(See also FigEJ) 



inverse of A, B = (A- 1 )' , which is given as 

f f f —f \ 
f f — f f 

f f Q f _f Q 

s -s s W 

-8 S 8 

s -s s ) 

with f — 1/r and s — 1/s. The generators of A* are given as the column vectors of B. We 
shall denote these generators as Qi (i = 1, 6), so that B = (Qi, Q2, Q3, Q4, Q5, Q&)- 

Since the TPMS's are projections of E(w) onto the 3D subspaces, E P , E D , and E G , 
their reciprocal lattices, A* p , A* D , and A^, are the intersections of A* along the relevant 3D 
subspaces. It can be shown that the sets of the three generators of these reciprocal lattices 



B = — 



1 
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FIG. 3: The orthogonal projections of the 6D lattice A for P- (left), G- (middle), and D-surfaces 
(right) are shown schematically in the complex plane of each component in z(w). The projection 
map of the lattice points of A are shown by open circles. The arrows indicate the directions of the 
orthogonal projections in each plane. 



are given by 

I Qi + Qs - Qg, Q2-Q4 + Qe, Q3 + Q4-Q5 } for A* P , 

I -Q2 + Q3 + <?4, Q1-Q3 + Qs, ~Qi + Q2 + Q 6 } for A* D , (8) 

I Qi - Qa - Qa, Q2-Q4- Q5, Q3 - Qb - Q& } for A* G - 

IV. ELECTRONIC BAND STRUCTURES OF THE TPMS'S 

For clarity, we shall focus on a specific situation in which the electron is constrained to 
move along the surface, which means that the electron cannot propagate through the region 
outside the surface. In order to realize the situation, an infinite potential well is introduced 
to confine the electron within a thin layer of constant thickness d built over the surface. 
Then by taking the limit of d — » + , one obtains the Schrodinger equation for the particle 
constrained to move along the surface . 3 i 4 i 5 i 6 Let us assume a curvilinear coordinates (g 1 , q 2 ) 
on the surface, for which the metric tensor = 1,2) is denned by ds 2 = J^-. g^-dgMg- 7 

where ds is an infinitesimal distance. The principal curvatures are denoted as k± — 1/R\ 
and k 2 = I/R2, with Ri and R 2 being the radii of curvature. The Schrodinger equation is 
written as, 

where g = det (#y) and (g^) = (gij)~ l . The first term describes the propagation of the wave 
function along the curved surface. The second term is the attractive potential energy asso- 
ciated with the local curvature,— resulting from the quantization procedure of the motion 
perpendicular to the surface. The normalization integral for ip is given by 

J J \iP(q\q 2 )\ 2 V9dqW, (10) 
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which is set to unity for square integrable wave functions. 

The above model, though seemingly too idealistic, is useful in answering fundamental 
questions about the quantum effects associated with geometrical surfaces. It may actually 
suit nano-structured materials composed of conductive thin curved films possibly to be 
realized in a near future. Since the model ignores the individual atomic potential energies, 
suitable systems should consist of substances with large Fermi wave length compared to the 
surface thickness or atomic distances, which may be achieved by semi-conducting substances 
(including graphene sheets). 

In order to study the electronic structures of the TPMS's, the WE representation must 
be incorporated to the above equation. Then one can study the relationship between the 
geometry, that is, the symmetry and topology, of the TPMS's and the electronic properties. 
This may be a useful approach to understand fundamental aspects of the problem for more 
general cases than minimal surfaces. Aoki et al. 10 formulated the Schrodinger equation with 
the stereographic projection, w = cot (9/2) exp (i(fi), into the following form: 

" W + cot ^ + W + m * ] = M 

Here the energy eigenvalue is represented by the dimensionless variable e = 8m£ 2 E/h 2 . e is 
a universal quantity among surfaces which just differ by scale, whereas the unit of energy, 
h 2 /(8m£ 2 ), depends on the scale (it is 9.525 meV for i = 1.0 nm, for instance). Remarkably, 
the above equation is free from the phase factor e ia (a is the Bonnet angle), so that it is 
identical for P-, G-, and D-surfaces. This is a result of the isometry of these surfaces as well 
as that the curvature potential depends only on the intrinsic (i.e., Gaussian) curvature for 
minimal surfaces, since — jjr-(/ci — K2) 2 = ^ K i K 2- Here we have used the vanishing of mean 
curvature, K\ + k 2 = 0. 

The topological differences among the surfaces should be taken into account through 
suitable choices of boundary conditions on the unit cell. Since the topology is crucial in re- 
stricting the propagation and interference of electronic wave functions, the electronic struc- 
tures should largely depend on it. It is an interesting question to elucidate how the band 
structures are affected by the topology. 

The Schrodinger equation can be solved numerically, by discretizing the coordinates 9 
and 4> and applying the finite difference method. In choosing the mesh values, special care 
is needed in the vicinity of the singularity at the branch points of the WE integrand to 
avoid the degradation of the resulting precision.— We employ a non-uniform mesh in which 
the maximal interval in the 3D embedding space would be minimized for a fixed number 
of mesh points. We calculate the band structures for P-, G-, and D-surfaces by using 4096 
mesh points per dodecagonal unit of Fig^ The results are shown in FigsJ3J|I3 and El 



V. INTERRELATION BETWEEN THE BAND STRUCTURES 

The crystallographic information given in Sec Mil provides a basis for understanding the 
interrelations between the electronic structures of these TPMS's. The boundary conditions 
for a dodecagonal patch of each surface is given as the phase factor multiplied upon trans- 
lating it to any of its neighbors. The phase factor is given by e 27rjk t for the translational 
vector t, where k is the 3D wave vector. Instead of the 3D vectors k and t, one can use the 
corresponding 6D vectors k and t which are restricted within the relevant 3D subspace. The 
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FIG. 4: The band structures of P-surface. The numbers shown on each of the energy levels indicate 
the corresponding IR's following the numbering scheme of Ref 27. The numbers are shown in black 
for special points and in gray for special lines. 
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FIG. 5: The band structures of D-surface. See the figure caption of Fig|I] 



phase factor is then rewritten as e 2mk4 , in which the vector t can safely be lifted to a 6D 
translational vector of H(u>) prior to the orthogonal projection. The boundary conditions 
can be generalized further by allowing the vector k, which should otherwise be restricted 
within the relevant 3D subspace for each surface, to take arbitrary values in the hypothetical 
6D "reciprocal space". Thus we are led to define the 6D band structures, e(k), which are 
periodic functions of the quasi-momentum, k = (ki, k 2 , k$, k^, fc 5 , fc 6 ). Since e(k) are derived 
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FIG. 6: The band structures of G-surface. See the figure caption of Fig^J 



from the periodicity of they are periodic functions of k with the periodicity of A*. It 

is understood that the band structures for each of the three TPMS's are the section of i(k) 
along the relevant 3D subspace. Hence, if kx is an arbitrary 6D vector along Ex (A = P,D, 
or G) and kx is the 3D counterpart, then the band structures for X-surface is given as 

It is readily shown that the intersection of a pair of subspaces, Ex and Ey (A, Y = 
P,D, or G), consists of just one point, which is at the origin. As a consequence, the 
energy eigenvalues of X-surface at T-point always coincide with those of Y-surface at T- 
point, because ex(0) = ey(0) = e(0). In this case, we can say that the band structures of 
X- and Y-surfaces "intersect" at T-point. 

A more general condition for the band intersection is given below. Recall e(k) has the 
periodicity of A*, so that 

£(k + J2n i Q i )=e(k), (12) 

i 

in which (i = 1, ...,6) are arbitrary integer coefficients. Therefore, if the two vectors kx 
in Ex and ky in Ey are connected through a scattering vector Yli n iQi lYl ^ '■> 

k Y = k x + y j n i Q h (13) 

i 

then the relevant eigenvalues coincide, i.e., ex(kx) = ey(ky). In particular, if Ylii n iQi e 
A* x © Ay, then one immediately finds that kx G A* x and ky e Ay. Since these wave vectors 
are at the center of the Brillouin zones (BZ's), the band intersection corresponds to T-point 
as discussed earlier. On the other hand, if ^ riiQi does not belong to A* x © A y , it is divided 
into two components in Ex and Ey which are not at the center of the BZ's. The latter case 
leads to band intersections at certain special points in the BZ's. The total number of band 
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P-surface 


D-surface 


G-surface 


r (t,h) 


r (t,r) 


r 


M (N 2 ) 


X (M, X) 




R (P 2 ) 




H 



TABLE I: The list of the special points for the band intersections, as derived in Appendix B 
The special points for the Bravais lattices of the dodecagonal unit patches are shown without 
parenthesis. The corresponding special points in the unfolded BZ's for P- and D-surfaces are 
shown by parenthesis. A superscript is added if two special points in the same star correspond to 
a single band intersection. 

intersections within the 6D BZ is given as the index // of the sublattice A^ © A Y in A*.— It 
can be shown that \x = 4 for the pair of P- and D-surfaces, 2 for that of P- and G-surfaces, 
and 1 for that of D- and G-surfaces. The corresponding band intersections occur at a set of 
special points in the BZ's of these TPMS's, as listed in TABLE U] (For the derivation, see 
Appendix B ) 

The coincidence of energy eigenvalues as the result of band intersections are observed in 
our numerical results as demonstrated in FigJ7| This phenomenon has been found previ- 
ously by Aoki et al . 10 ' 11 in an empirical way. However, our theory provides a more clear 
account for this phenomenon in terms of the 6D algebra of Bonnet transformations. It also 
gives the exhaustive list of possible band intersections in TABLE [I] Moreover, the idea of 
band intersections can be applied generally to the band structures of TPMS's connected via 
Bonnet transformations. 

VI. SYMMETRY ANALYSIS OF THE BAND STRUCTURES 

Symmetry has various implications in understanding the properties of materials involving 
the quantum mechanics of electrons. Dynamically, it implies the selection rules for quantum 
transitions or various forms of conservation laws. While statically, symmetry serves as 
the classification principle of structures and also provides basic clues for understanding the 
formation (or stability) of a particular material. It should also have significant implications in 
understanding the detailed characteristics of the present systems. Therefore, in the following, 
we shall carry out a detailed analysis on the symmetry aspects of the electronic structures 
of the TPMS's. 

A subspace of Bloch states 0k at a wave vector, k, with an identical energy eigenvalue 
can be assigned an irreducible representation (IR) of G^, i.e., the group of k. In practice, 
the representation matrices are obtained numerically, and are then analyzed by the use 
of the orthogonal relationships of IR's . 25 i 26 In fact, the discretization of the surfaces used 
for the numerical analysis 10 diminishes the space group, G, to P4/mmm, P^/nnm, and 
I4i/acd for P-, D-, and G-surfaces, respectively. Hence, for each surface, we can numerically 
analyze the IR's only at k-points whose group G^ are subgroups of the diminished space 
group. In addition, we may use the band branching rules (for instance, T 9 always splits into 
A3 and A5 for the space group Im3m) to determine the IR's for the rest of the k space. 
The numbers indicated on the band structure functions in FigsJU and El represent the 
relevant IR's of G^; we comply with the numbering scheme given in Refl27l These IR's 
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FIG. 7: Energy diagrams showing the coincidence of energy levels at the special points listed in 
TABLE QJ The superscript for each special point symbol (e.g., N 6 ) indicates the number of the 
equivalent points in the BZ, constituting the star. Each energy level is assigned the relevant IR, 
whose dimensionality is shown by parenthesis. 



describe the symmetry features of the relevant eigenstates. They also serve as criteria for 
the degeneracies in the band structures. This is useful since otherwise degenerate bands can 
slightly fall apart as artifacts of the diminishing of symmetry involved by the discretization 
procedure. 



A. Band sticking 

Degeneracies occur due to several symmetry reasons. At high-symmetry special points 
in the BZ, we often observe several band energy curves sticking together into one energy 
level. The degeneracy of such an energy level is given by the dimensionality of the relevant 
IR. For the case of non-symmorphic space groups, the degeneracies at a special point on the 
boundary of the BZ can be larger than symmorphic counterparts when the transformation of 
eigenstates by is represented as a multiplier group, 25 In the band structures of G-surface, 
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for instance, we may observe relatively large degeneracies at H- and P-points; they are 6- 
fold (Hi), 4-fold (Pi), and 2-fold (H2, H3, H4, P2, P3). Furthermore, we may find additional 
degeneracies due to the time-reversal symmetry on these special points . 25 i 28 Namely, an 
energy level of the IR H 3 always sticks with one of the IR H 4 , whereas the same holds for 
the pair of the IR's, P 2 and P 3 , both resulting in 4- fold degeneracies. Note the space group 
for G-surface, Ia3d, is highly non-symmorphic including four-fold screw axes and two types 
of glide planes (one type is d-glide). 



B. Nodal lines 

An interesting feature of the present systems is that, due to the 2D character of the 
surfaces, some of the eigenstates have nodal lines, i.e., lines along which an eigenstate has zero 
amplitude. Here, we shall study the occurrence of nodal lines in somewhat detail from the 
symmetry standpoint, because nodal lines have the following significant implications for the 
electronic properties. Firstly, the density of nodal lines is correlated with the magnitude of 
energy eigenvalue. Secondly, nodal lines occur as a consequence of the quantum interference 
of an electronic wave and are closely related to the symmetry and topology of the surfaces. 

There is an important class of nodal lines which can be understood in terms of the 
symmetry properties of the eigenstates; a nodal line should appear when the relevant IR 
has odd parity with respect to a symmetry element which leaves the positions of the nodal 
line invariant. There are two types of such symmetry elements: mirrors (a nodal line will 
appear at the intersection between the surface and a mirror plane) and two-fold axes lying 
in the surface (a nodal line will appear along such an axis). 

In the space group Im3m of P-surface, the symmetry elements related to the generation 
of nodal lines are 

• (two-fold axes U d ) U xy , U Sy , U yz , U yz , U zx , U zx , 

• (mirrors ad) a xy , a xy , a yz , a yz , a zx , a zx , 

• (mirrors a) a x , a y , a z . 

Note that the elements of the factor group are implied by the symbols of point group 
elements. Throughout this paper, the symbols of point group elements comply with those 



used in Ref|27i One can find that the above three types of symmetry operations leave 
the high-symmetry geodesies on P-surface, i.e., a-, j3-, and 7-type geodesies as defined in 
Fig |2J invariant. The correspondence between the types of symmetry operations and those 
of geodesies are given by a «-> Ud, (3 <-> ad, and 7 <-> a (see FigHJ). 
A similar list for D-surface (space group Pn3m) is given as 

• (mirrors ad) a xy , a xy , a yz , a yz , a zx , a zx , 

• (two-fold axes U d ) U xy , U xy , U yz , U yz , U zx , U zx , 

• (two-fold axes U) U x , U y ,U z . 

In this case the correspondence is established as, a «-> a d , (3 <-> Ud, and 7 «-> U (see FigJTJ). 

On the other hand, no such symmetry elements exist for G-surface, that is, there are no 
mirror planes (though there are several glide planes) nor two-fold axes lying in the surface. 
Hence no nodal lines are predicted by symmetry for G-surface. 
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k-point 


G k /T 


IR (nodal lines) 


r (0,0,0) 


o h 


r 2 (6C/ rf , 6a d ), r 6 (3cr, 6cr d ), T 7 (6U d , 3a),T 8 (3a) 


A (^,0,0) 


G± v 


A 2 (2a, 2o- d ),A3(2a d ),A4(2(7) 


H (1,0,0) 


o h 


iJi(6C/ d ) ; H 2 (6a d ), H 6 {6U d , 3a, 6a d ), H 7 (3a), H 8 (3a) 


S (k x , k x , 0) 




^2(U d , a d ), S 3 (cr, cr d ), E 4 (C/ d , a) 


G {k x , 1 — k x , 0) 


Civ 


Gi(U d ),G 2 {a d ), G 3 (U d , a, a d ), G A {a) 


iV (±,±,0) 


D 2h 


N x {u s y), N 2 (a s y, a z ), N 3 (u x y , u*y, a z , a x y), N 4 (u x y, a x y, a s y) 
N 5 (u x y,a z ,a x y,a x y),N 6 (a x y),N 7 (u x y,u x y,a x y),N 8 (u x y,a z ) 


A (fej;, fej;, fej;) 


Czv 


A 2 (3a d ) 


F (fc x , 1 - k x , 1 - k x ) 


C^v 


F 2 (3a d ) 


D (5, 2: &z) 


Civ 


D 2 (a x y),D 3 (a x y,a x y),D 4 (a x y) 


P (- - -) 

1 \ 2 ' 2 ' 2 / 


T d 




S (fej; , fcy, 0) 


C s 


S 2 (a) 




C s 


e 2 (a d ) 


(A^, 1 — k x ,k z ) 


C s 


S 2 (cT d ) 



TABLE II: The list of the IR's for P-surface (Im3m) entailing nodal lines by symmetry. The IR's 
shown have odd parity with respect to the symmetry elements given by parenthesis; these elements 
are of the types, U d , a d , and a, which are responsible for the generation of nodal lines along a-, 
/?-, and 7-geodesics, respectively. A set of symmetrically equivalent elements are shown by the 
numbered type symbols; e.g., 3a = {a x ,a y ,a y }. 



It is deducible from the character table of a given k-point in the BZ whether a given IR 
entails any nodal lines by symmetry. One can consult the character tables for the space 



groups in Refl27l It is necessary from the previous paragraph that the group of k, G^, 
contains at least one symmetry element (s) responsible for the generation of nodal lines. 
Then, if an IR of has odd parity with respect to those symmetry elements, a nodal line 
will appear accordingly. One should bear in mind that the representative of the relevant 
factor group (7(c) (7 is the point group operation and c the additional translation, which 
is usually taken to be zero when the space group is symmorphic) does not necessarily leave 
the position of the nodal line invariant. Instead, one should consider (/y|r) with 

r = q-7q = c-R4 7|c) , (14) 

where q is any position in the relevant mirror or axis. Here, Rf |c) = ( 7 |c) q — q is always a 
Bravais vector. The relevant matrix for (7|r) is given by 

D[(7|r)] = cxpHkR^] J D[( 7 |c)]. (15) 

Therefore, in order to discern the parity of an IR at the position of a nodal line, the char- 
acter for the representative element of the relevant factor group (the value given in the 
character table) should be multiplied by exp [— ikR^ 7 ' - 1 ]. The exhaustive lists of IR's which 
entail nodal lines by the above argument are given in TABLES |H] and IIHI for P- and D- 
surfaces, respectively. For each of the IR's listed, the corresponding nodal lines are given by 
parenthesis. 
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k-point Gk/T 



IR (nodal lines) 





M 1 (U z ),M 2 (U z ,2a d ), M 4 (U x y, [/*») 



r 2 {QU d ,6a d ),T 6 (6a d ),T 7 {6U d ) 
A 2 (2a d ),A 3 (2a d ),A 5 (U) 

x 2 (2a d ),x 3 (u*, xjy\ uy z ),x 4 (u x ) 

£2(0^, U d ), E 3 (a d ),T, i (U d ) 



&2(a d ) 



TABLE III: The list of the IR's for D-surface (Pn3m) entailing nodal lines by symmetry. The IR's 
shown have odd parity with respect to the symmetry elements of the types, a d , U d , and U, which 
are responsible for the generation of nodal lines along a-, (3-, and 7-geodesics, respectively. 

For G-surface, we may still find related nodal lines for eigenstates on the special points 
T and H, corresponding to band intersections. These eigenstates are in one-to-one corre- 
spondence with eigenstates of P- and D-surfaces on the relevant special points (see also a 
discussion in Sec lVII)) . and they are identical intrinsically within the dodecagonal unit patch. 
Therefore if some nodal lines are present for symmetry reasons for P- and D-surfaces, the 
corresponding states of G-surface should also have nodal lines. Illustrative examples are 
given in Fig|HJ which shows two sets of eigenstates with IR's {H 7 , R 7 , T 2 } and {H 2 , T 7 , T 7 } 
for {P, D, G}-surfaces, respectively. One can see that the sets of nodal lines are intrinsically 
the same within each set. 

We find that some eigenstates have nodal lines which do not correspond to any of the 
three types of high symmetry geodesies. This kind of nodal lines are by no means easy to 
predict because they are not directly related to mirrors nor two-fold axes in the surface. 
However, such nodal lines would appear only at sufficiently high energies. 

VII. DISCUSSIONS 

The idea of band intersections discussed in SecjV] can be paraphrased in the following 
way. The Bloch condition for the eigenstates of a crystal is given as the proper phase 
relationships of the reference unit cell to its neighbors. For the present case, we have shown 
that the phase relationships of a dodecagonal unit to its twelve neighbors become identical 
between two different surfaces at the k points of a band intersection. Therefore, the partial 
Hamiltonians acting on the subspaces of states at these k points will become equivalent. 
The coincidence of energy eigenvalues is just a result of this fact. If there are several members 
in the star of the special points, all these members give a set of symmetrically equivalent 
phase relationships of the dodecagonal unit patch; for instance, the stars of M-point (s.c.) of 
P-surface and X-point (f.c.c.) of D-surface forming band intersections are three membered. 
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P-surface 




H-point 3rd energy level (IR-H ? ) H-point 5th energy level (IR-H 2 ) 



□-surface 




R-point 3rd energy level (IR-R 7 ) f-point 6th energy level (IR-r,) 




r-point 6th energy level (IR-ry r-point 10th energy level (IR-r ? ) 

FIG. 8: (Color) In each column eigenstates with identical energy eigenvalue corresponding to a 
band intersection are shown. The dots are the actual mesh points for the numerical calculation. 
Opposite signs of the wave functions are indicated by the blue and red colors, while the black dots 
indicate the amplitude whose absolute value is close to zero, corresponding to the nodal lines. The 
states in the set {(a), (b), (c)} have nodal lines along all the 7 type geodesies, while ones in the set 
{(d), (e), (f)} have nodal lines along all the (3 type geodesies. The 7 type geodesies for G-surface 
in (c) constitute the helical lines of both chirality along x, y, and z directions . 



Except for at the band intersections, the phase relationships cannot coincide among the 
three surfaces because of their topological differences. 

The equivalence of partial Hamiltonians between two surfaces at a band intersection 
implies a one-to-one correspondence between the subspaces of states of the same energy 
level. If the two sides of the surfaces are distinct, the correspondence would be simple so 
that one could expect a one-to-one correspondence between the IR's between the surfaces. 
However, as we see in FigJTl the correspondence is somewhat obscure due to the reduction 
of the primitive unit cells for P- and D-surfaces when the two sides are indistinguishable. 
The unfolding of the BZ's for these two surfaces results in the splitting of a single band 
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intersection in the 6D BZ into two different special points (see TABLE U}. The resulting two 
special points can either be (i) equivalent, as in the case of two iV-points for P-surface, or (ii) 
inequivalent, as in the case of T- and if -points for P-surface. For the case (i), a subspace of 
states at a band intersection is divided into two subspaces associated with the two k points 
in the same star, so that the dimensionality of the IR's is reduced by half from that of the 
original subspace. Meanwhile, in the case (ii), a subspace of states at a band intersection is 
mapped onto a single subspace associated with either of the different special points. This 
argument is based on that each subspace at a single k point is in general associated with a 
single IR, and that the total dimensionality of the original subspaces at a band intersection 
is identical between the relevant two surfaces. 

We expect that there is a certain connection between the IR's of different surfaces at their 
band intersections. It appears that an IR of P-surface corresponds uniquely to an IR of D- 
surface, while an IR of G-surface corresponds to two IR's of P-surface as well as D-surface. 
In FigJ7| one finds that the IR r 9 of G-surface corresponds to either r 9 or P 4 of P-surface 
as well as to either T 9 or P 4 of D-surface. Interestingly, the connection between the IR's of 
P-surface at P-point and those of G-surface at P-point can be uniquely determined by the 
dimensionality of the IR's as 

Pi(l),P 2 (l) <- # 2 (2), 

P 3 (2) <-> P 3 + P 4 (2 + 2), 
P 4 (3),P 5 (3) <-> Pi (6), 

where the dimensionality of each IR is denoted by parenthesis. (Recall that the IR's P3 and 
P 4 for G-surface is degenerate due to the time-reversal symmetry.) For rigorous proofs of 
the general relationships between the IR's of the three surfaces, more detailed investigations 
on the group theoretical aspects of the problem may be required. Still, our results can be a 
starting point for further studies. 

The relationship between the symmetry properties of eigenstates and the occurrence of 
nodal lines are established in Sec JVIl The existence of nodal lines in an eigenstate necessitate 
a sufficient amount of kinetic energy, so that IR's with a large number of nodal lines cannot 
appear in a lowest part of the energy spectrum. In TABLES ITT1 and ITTT1 one can find that the 
following IR's accompany a significant number of nodal lines: {T2, Tq, T?, Hq} of P-surface 
and {T 2 , Pi, P2, Rq} of D-surface. These IR's are naturally absent in FigsJU and 03 It is 
also seen that IR's with nodal lines appear roughly in the increasing order in the number 
of nodal lines as one goes up the energy axis. Therefore, nodal lines can be regarded as 
an important feature of eigenstates present in real-space in order to understand the overall 
characteristics of the systems. 

For a full account of the electronic structures of the present systems, a more detailed 
study on the topological aspects of the problem should be important. Due to the topolog- 
ical characters of the surfaces, there may be a "topological invariant" associated with the 
continuum of eigenstates along a dispersion curve in the band structures. Such a topological 
invariant might be given as the phase winding numbers associated with characteristic closed 
paths on the surface. The topological characters of the eigenstates may explain the connec- 
tivity aspects of the energy bands. Further properties of the energy bands (e.g., the flatness, 
particle-/hole-type dispersions) might also be explained based on real-space pictures of the 
relevant eigenstates. 

It is instructive to note that the arguments throughout this paper are not binded by the 
particular form of the curvature potential as long as it is determined intrinsically. Hence, by 
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generalizing the potential to any function of the Gaussian curvature one may be able to deal 
with secondary effects of curvatures, such as the effect of varying thickness of the surface 
depending on the curvatures. 

In summary, we have studied in detail the band structures of the strongly constrained 
electron systems onto P-, D-, and G-surfaces. The crystallographic relations of these surfaces 
in six dimensions have proved to be useful in understanding the interrelations between the 
electronic structures of these surfaces. By introducing the hypothetical 6D band structures, 
the coincidence of eigenstates among different surfaces at particular sets of k points has been 
explained in terms of the idea of band intersections. We have given the exhaustive list of 
band intersections for the three surfaces. A numerical investigation of the band structures 
with a detailed analysis of the symmetry properties has also been performed. We have 
assigned IR's of the space groups to the energy bands displayed in Figs J4l5l and El As a key 
feature to understand the consequences of the symmetry of the eigenstates, the occurrence 
of characteristic nodal lines have been discussed somewhat in detail. We have given the 
conditions for nodal lines to be predicted by symmetry and also the exhaustive lists of IR's 
which entail nodal lines for P- and D-surfaces. 
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APPENDIX: THE INDEX OF A y IN A* 

We consider the sublattices of A* given in the form © Ay (X, Y = P, D, or G). From 
Eq.fjHJ), it is straightforward to obtain the following: 

Ap © A* D = | i (4 + 4 - Q 6 ) + j (4 - 4 + 4) + & (4 + 4 - 4) 

+ I (-4 + 4 + 4) + ™ (Q a - Q 3 + 4) + n (-4 + 4 + 4) } ,(A.l) 

Ap © Aq = { % (Qi + 4 - 4) + j (4 - + Qe) + k (Q 3 + Q*~ Qs) 

+ l{Qi-Qi- Qe) + m (Q 2 - g 4 - 4) + n (4 - 4 - 4) } , (A.2) 

a* d ®a* g = { % (-4 + 4 + 4) + j (4 - 4 + 4) + k (-4 + 4 + 4) 

+ 1 (4 - 4 - 4) + m (4 - 4 - 4) + n (4 - 4 - 4) } , (a.3) 

where i, j, k, I, m, and n are arbitrary integers. In the following, the orthonormal basis sets 
of the 3D subspaces Ep, Ed, and Eg are denoted as ep t i, eo,i and ec,i (i = 1,2,3), respec- 
tively. These basis vectors are obtained as the row vectors of the corresponding orthogonal 
projectors, lip, Up, and 11^. 
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i. a* p © a* l 



Ap © A* D is the sub lattice of A* = {Y^a = i n iQi} for which the coefficients are given by 



/nA 




( 


i 






( 


1 











1 


-1 


\ 








J 











1 





-1 





1 






= M PD 




k 
I 




, M PD = 










-1 


1 
1 


1 
1 


-1 










n 5 






in 








1 





-1 





1 







\n 6 / 




\ 


n 


) 




\ 


-1 


1 











1 


/ 



(A.4) 



Since the index | detMp£>| = 4, A* is divided into four sublattices equivalent to A* P © A* D . 
Each sublattice is specified by a pair of indices, (p, q): 

p = rii + n 3 + 77-4 + tiq mod 2, (A. 5) 

q = n 2 + n 3 + n 5 + n 6 mod 2. (A. 6) 

In particular, the sublattice A* p Q)A* D is specified by {p = 0, q = 0). It gives rise to the band 
intersection at T-point, kp = kp> = 0, in Eq. (fH?|) . The other three sublattices are obtained 
by translating A* p © A^ by Qi (p = 1, q = 0), Q 2 (p = 0, q = 1), and Q 3 (p = 1, q = 1). For 
each of these sublattices, the vectors kp and kp, in Eq. (jl3j) can be taken as 

k P = (e Pi2 + e Pi3 )/(2r£), k D = (e D , 2 - e Dt3 )/{2s£) 

forp=l,g = (k P -k D = Qi), (A.7) 

fcp = (e P ,i + e Pt3 )/(2r£), k D = (-e A1 + e D>3 )/(2s£) 

{orp = 0,q = l (k P -k D = Q 2 ), (A.8) 

fcp = (e P ,i + e P , 2 )/(2r£), fc D = (e D) i - e Dj2 )/(2s£) 

forp=l,g=l (fc P -fc D = Q 3 ). (A.9) 

It is readily shown that these three band intersections are observed at the three membered 
star of M-point in the BZ of P-surface as well as at the three membered star of X-point in 
the BZ of D-surface, where the dodecagonal patches are taken as the primitive cells. 



2. A* P © A* G 



A similar argument holds for Ap © A^, for which 

f 1 1 \_ 
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1 1 
0-11-1-10 
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V — i i 0-10-1/ 



M 



PG 



(A.10) 



The index is |detMp G | = 2. Hence, there are two equivalent sublattices specified by an 
index, r: 



n 4 + n 5 + n e mod 2. 



(A.ll) 
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Apart from the trivial band intersection for r = at r-point, we have another band intersec- 
tion for the sublattice r — 1, which is obtained by translating A* P © A G by Q4. The relevant 
vectors kp and k G can be taken as 

k P = (e P)1 - e P)2 + e Pi3 )/(2r£), £ G = e G ,i/(2M) 

forr = l (kp — kc = Qa)- (A.12) 

The band intersection is therefore observed at R-point in the BZ of P-surface as well as at 
H-point in that of the G-surface, where the dodecagonal patches are taken as the primitive 
cells. 



A* D e A* c 
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The index is | det M^g | = 1 . Therefore A 
other than the trivial one at r-point. 



D 



)A G 



A*, and there is no extra band intersection 
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